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Abstract 

The fortran version of the AcerDET package has been published in [1], and used in the 
multiple publications on the predictions for physics at LHC. The package provides, starting 
from list of particles in the event, the list of reconstructed jets, isolated electrons, muons, 
photons and reconstructed missing transverse energy. The AcerDET represents a simplified 
version of the package called ATLFAST, used since several years within ATLAS Collabora¬ 
tion. In the fast simulation implemented in AcerDET, some functionalities of ATLFAST are 
absent, but the most crucial detector effects are implemented and the parametrisations 
are largely simplified. Therefore it is not representing details neither of ATLAS nor CMS 
detectors. 

This short paper documents a new C++ implementation of the same algorithms as 
used in [1], We believe that the package can be well adequate for some feasibility studies 
of the high px physics at LHC and at planned ppFCC. The further evolution of this code 
is planned. 


[1] E. Richter-Was, AcerDET: A Particle level fast simulation and reconstruction pack¬ 
age for phenomenological studies on high p(T) physics at LHC , hep-ph/0207355. 


1 Package is available for the web page http://erichter.home.cern.ch/erichter/AcerDET.html 



PROGRAM SUMMARY 


Title of the program: AcerDET version 2.0 
Operating system: Linux 
Programming language: C++ . 

External libraries: HepMC, ROOT. 

Size of the compressed distribution directory: about 25 kB. 

Key words: Fast simulation, Physics at LHC. 

Nature of physical problem: A particle level fast simulation and reconstruction package. 
The package provides, starting from the list of particles in the event, list of reconstructed 
jets, isolated electrons, muons and photons and reconstructed missing transverse energy. 
The algorithmic part is just rewritten from fortran version of the package AcerDET-1.0 

[1], The default input is generated event in format of HepMC structure [2], output is 
in form of ROOT file [4], containing tree with reconstructed objects and set of control 
histograms. Distribution version includes also example of the main program for execution 
with PYTHIA 8.2 generator [3]. Implemented set of parametrisations is not representing in 
details ATLAS or CMS detectors, some of them are simple and can be considered rather 
as place-holders for future adaptation of any detector. Nevertheless, we believe that the 
package will be well adequate for some feasibility studies on the high px physics at LHC. 

[1] . E. Richter-Was, AcerDET: A Particle level fast simulation and reconstruction 
package for phenomenological studies on high p(T) physics at LHC , hep-ph/0207355; 
http://erichter.web.cern.ch/erichter/AcerDET.html 

[2] . M. Dobbs and J.B. Hansen, HepMC a C++ Event Record for Monte Carlo Gen¬ 
erators. , Comput. Phys. Commun. 134 (2001) 41; 

http://lcgapp.cern.ch/proj ect/simu/HepMC/ 

[3] . T. Sjstrand, S. Mrenna and P. Skands, A Brief Introduction to PYTHIA 8.1 
Comput. Phys. Comm. 178 (2008) 852; T. Sjstrand, S. Mrenna and P. Skands, PYTHIA 
6.4 Physics and Manual, JHEP05 (2006) 026; 

http://home.thep.lu.se/ torbjorn/pythia81html/Welcome.html 

[4] . ROOT: Data analysis framework 
https://root. cern. ch/drupal/ 
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1 Introduction 


The potential of the LHC detectors for physics at high px is very rich, as first three years 
of data taking during so called Run I have shown. 

The discovery of the Higgs boson, searches results for Supersymmetry particles, Exotic 
particles, New Vector Bosons are very impressive. All these is thanks to the high sensitivity 
of the detectors in terms of acceptance and identifying efficiencies to variety of signatures: 
photons, electrons, muons, multi-b-jets, r-jets, missing transverse energy. 

Those high sensitivities allow to study very exclusive signatures; the discovery potential 
is limited in some cases by the rare background processes. It is evident, that the multi¬ 
jet and multi-b-jet production in association with known vector bosons, W and/or Z, 
or with top-quark pair is the the serious background to the several observables as well. 
The importance of effect beyond the leading order: effects coming from the finite width 
or angular spin correlations (eg. from the intermediate resonances decays) is becoming 
more evident. Limitations and complementarity of the matrix element and parton shower 
predictions for the diversity of signatures need to be investigated. 

The package for particle-level simulation and reconstruction is one of the intermediate 
steps between simple parton-level analysis and very sophisticated and CPU consuming 
full detector simulation. The package provides, starting from list of particles in the event, 
list of reconstructed jets, isolated electrons, muons and photons and reconstructed missing 
transverse energy. It can serve for phenomenological feasibility studies on the prospects 
for observability of a given signature. Such package can be useful for the dedicated com¬ 
parisons between matrix element and parton shower predictions. In that case, one can 
compare experimental signatures (reconstructed jets, leptons, photons) to quantify size of 
the discrepancies between different predictions in a straightforward way. 

The Fortran version of presented here package [1] has played such role since several 
years, more than 80 papers published in the last 10 years was using it or this purpose [2], 
Here, we present the new C++ implementation, where all the algorithmic components have 
been preserved, changed is format of the input event. Now comonly accepted HepMC 
and the output format of ROOT tree with reconstructed objets and control histograms is 
introduced. 

The package simulates some key features of the LHC detectors like ATLAS and CMS. 
It is based on the calorimetric energy deposition for jets reconstruction and tracks re¬ 
construction for electrons and muons. It takes into account very high granularity of the 
electromagnetic calorimetry for the photon reconstruction. The missing transverse energy 
is calculated for the total energy balance of the reconstructed objects. The capability for 
the identification of b-jets and r-jets is also explored. Implemented set of parametrisations 
is not representing in details performance of neither ATLAS nor CMS detectors. Never¬ 
theless, we believe that the package is still adequate for feasibility studies on the high px 
physics at LHC and offer starting option also for future FCC program. 

The paper is organised as follows. In Section 2 we discuss (recall) algorithms used 
for objects reconstruction and show benchmarking distributions. Given the Higgs boson 
discovery at mass 125 GeV, benchmark which were published in [1] are now updated to 
this value of teh Higgs mass. Also Pythia 8.1 showering model used here are different 
than what was used for [1], hence some benchmark figures are different despite same 
algorithms used for fast simulation and reconstruction. Section 3 gives an outlook. We 
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don’t recall comparison with ATLAS performance numbers (Appendix A of [1]) as those 
are not updated as well with the better inderstanding of the detector using Run I data 
and are much more now detailed and analysis dependent. One should just to look for the 
publications related to the signature of interest. It was however checked that numbers 
from Appendix A of [1], are reproduced by AcerDET-2.0 code. 

In the Appendixes A-F, we give more technical details concerning input parameters 
and output structure. In Appendix G we show control output logfile. 

2 Simulation and reconstruction 

Fully or partially generated event i.e. event generated including or not QED/QCD initial 
and final state radiation, fragmentation, hadronisation and decays of unstable particles 
can by analysed by this package. The list of the partons/particles of the generated event 
should be stored in the HepMC format. It is then copied into internally used vector of 
particles, with added few flags i.e. additional information regarding particles status and 
its type. This information is used by the AcerDET algorithms for event-be-event simulation 
and reconstruction. 

Events simulation is limited to the following steps: 

• Deposition of the particles energies in the calorimetric cells. 

• Smearing of the energy of electrons, photons, muons with the parametrised resolu¬ 
tions. 

• Smearing of the energy of hadronic clusters and not-clustered cells with the parametrised 
resolution. 

Events reconstruction is limited to the following steps: 

• Reconstruction of the calorimetric clusters. 

• Verification of the isolation criteria for electrons/photons/muons. 

• Rejection of clusters associated with electrons and photons. 

• Acceptation of remaining clusters as hadronic jets and identification (labeling) of 
those associated with b-quarks, c-quarks, tau-leptons. 

• Jets energy calibration. 

There is no clear separation between the simulation and reconstruction parts in the struc¬ 
ture of algorithms. A larger blocks (algorithms) may realise both tasks. 

As a final result provided are four-momenta of reconstructed objects: electrons, pho¬ 
tons, muons, labeled and calibrated jets and calculated is total missing transverse energy 
in analysed event. 
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2.1 Calorimetric clusters 


The transverse energy of all undecayed particles stored in the event record, 

(InputRecor: Particles, except of muons, neutrinos and other invisible particles 2 eg. 
the lightest SUSY particle, are summed up in the map of calorimetric cells with a given 
granularity in (?y x cf) coordinates (default: 0.1 x 0.1 for |^| < 3.2 and 0.2 x 0.2 for |^| > 3.2, 
with the calorimetric coverage up to |?y| = 5.0). As an effect of the solenoidal magnetic 
field in the inner part of the detector we assume that the 4 > position of charged particles 
with transverse momenta above the threshold (default: px > 0.5 GeV) will be shifted as 
parametrised in function Particle: : f oldPhi. The contribution from all charged particles 
with transverse momenta below that threshold is neglected. 

All calorimetric cells with the transverse energy greater than a given threshold (de¬ 
fault: Et > 1.5 GeV) are taken as possible initiators of clusters. They are scanned in 
order of decreasing Ex to verify whether the total Ex summed over all cells in a cone 
A R = \J A 2 ?y + A 2 cj) exceeds the minimum required threshold for the reconstructed clus¬ 
ter (default Ex > 5 GeV). Cells with deposited transverse energy below the threshold 
(default: Ex = 0) are not accounted for. As a coordinates ( rj clu x cj) clu ) of the recon¬ 
structed cluster taken are the coordinates of the bary-center of the cone weighted by the 
cells Ex for all cells inside the cone around the initiator cell. Information about energy 
deposition in cells is stored in OutputRecord: Cells 

All reconstructed clusters are stored in the container of ClusterData objects, 
OutputRecord:Clusters. Fig. 1 shows the A p (left) and Af> (right) distribution between 
the reconstructed bary-center of particles falling within the geometrical cluster cone and 
the reconstructed cluster position. 




Figure 1: The A p (left) and A <f> (right) distribution between the reconstructed bary-center 
of particles falling within the cluster cone and the reconstructed cluster position. Shown 
are results for generated gg —>• H, H —>• uu process with ran = 125 Ge V. 

2 User may want to make any particle invisible to the detector. For this one should overwrite its pdgld 
code in the vector<Particle> to that specified for the invisible particles in acerdet.dat file. 


4 




2.2 Isolated muons 

Algorithm reconstructing isolated muons uses information of the generated muons, recon¬ 
structed calorimetric clusters and the cells map. 

Isolated muon candidates are searched for in the container InputRecord.particles. 

The inverse muon four-momentum is smeared according to the Gaussian resolution parametrised 
with function Smearing: :forMuon (default: a = 0.05% -pr)- The muon direction remains 
unsmeared. 

For all muons which pass selection criteria in px and r ] (default: px > 6 GeV and 
|?y| < 2.5 GeV), isolation criteria, in terms of the distance from calorimetric clusters 
and of maximum transverse energy deposition in cells in a cone around the muon, are 
then applied (defaults: separation by A R > 0.4 from other clusters and &T < 10 
GeV in a cone A R = 0.2 around the muon). All muons passing the isolation criteria 
are stored in the container OutputRecord.Muons and those not passing are stored in 
OutputRecord.NonisolatedMuons. 

As a control physics process the gg —>• H —> ZZ* —> production with the Higgs 

boson mass of 125 GeV is used. The default isolation selection has 97.8% efficiency for 
muons passing kinematical selection with pjt ’ /12 > 20 GeV and pj?’^ 4 > 7 GeV. As the 
predicted intrinsic width of the Higgs boson of this mass is very small, the resolution is 
completely dominated by the resolution assumed for the single muon reconstruction. Fig. 2 
(left) shows the reconstructed distribution of the 4-muon system. The assumed single muon 
transverse momenta resolution of a = 0.05% -px leads to the a m = 1.57 GeV resolution for 
the invariant mass of the four-muon system originating from the H ^ ZZ* ^ Ap decay. 



Figure 2: Left: The multiplicity of hard-process isolated muons (left) and reconstructed 
isolated muons (right) for gg —>• H, H —>• ZZ* —>• Ap events with mH = 125 GeV. 

2.3 Isolated electrons 

Algorithm reconstructing isolated electrons uses information of the generated electrons, 
reconstructed calorimetric clusters and the cells map. 

Isolated electron candidates are searched for in the InputRecord.particles. The 
electron four-momentum is smeared according to the Gaussian resolution parametrised 
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with function Smearing: :forElectron (default: a = \2%/\fE). The electron direction 
remains unsmeared. 

For all electrons which pass selection criteria in px and g (default: Pt > 5 GeV and 
|?y| < 2.5 GeV), the associated reconstructed calorimeter cluster is identified (default: 
AR ej cluster < 0.1). Electron isolation criteria, in terms of the distance from other clusters 
and of maximum transverse energy deposition in cells in a cone around the electron, are 
then applied (defaults: separation by A R > 0.4 from other clusters and Y) Rt < 10 GeV 
in a cone A R = 0.2 around the electron). All electrons passing the isolation criteria are 
stored in OutputRecord. Electrons and the clusters associated with them are removed 
from the OutputRecord.Clusters. 

As a control physics process the pp —> —>■ ZZ* —* 4e production with the Higgs 
boson mass of 125 GeV is used. The default isolation selection has 95.2% efficiency for 
electrons passing kinematical selection with pjl ’ 62 > 20 GeV and p^' Li > 7 GeV. As the 
predicted intrinsic width of the Higgs boson of this mass is very small, the resolution is 
completely dominated by the resolution assumed for the single electron reconstruction. 
Fig. 3 (left) shows the reconstructed distribution of the 4-electron system. The assumed 
single electron energy resolution of a = 12%/\/E leads to the a m = 1.26 GeV resolution 
for the reconstructed invariant mass of the four-electron system originating from the H —>• 
ZZ* —>• 4e decay. 



Figure 3: Left: The multiplicity of hard-process isolated electrons (left) and reconstructed 
isolated electrons (right) for gg —>• H, H — >• ZZ* —>• 4e events with m# = 125 GeV. 

2.4 Isolated photons 

Algorithm reconstructing isolated photons uses information of the generated photons, 
reconstructed calorimetric clusters and the cells map. 

Isolated photon candidates are searched for in the InputRecord.particles. The 
photon four-momentum is smeared according to the Gaussian resolution parametrised 
with function Smearing: :forPhoton (default: a = 10%/V^E). The photon direction 
remains unsmeared. 

For all photons which pass selection criteria in px and r] (default: px > 5 GeV and 
|?y| < 2.5 GeV), the associated reconstructed calorimeter cluster is identified (default: 
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A R-y,cluster < 0.1). Photon isolation criteria, in terms of the distance from other clusters 
and of maximum transverse energy deposition in cells in a cone around the photon, are 
then applied (defaults: separation by A R > 0.4 from other clusters and Y) Et < 10 GeV 
in a cone A R = 0.2 around the photon). All photons passing the isolation criteria are 
stored in OutputRecord.Photons and the clusters associated with them are removed from 
the OutputRecord.Clusters. 

As a control physics process the Standard Model gg —>• H —>• 77 production with the 
Higgs boson mass of 125 GeV is used. The isolation selection has 98.0% efficiency for 
photons passing the kinematical selection of gift > 40 GeV and p j? > 25 GeV. As the 
predicted intrinsic width of the Higgs boson of this mass is much below 1 GeV, the resolu¬ 
tion of the reconstructed invariant mass of the di-photon system is completely dominated 
by the resolution assumed for the single photon reconstruction. Fig. 4 (left) shows the 
reconstructed distribution of the di-photon system of photons passing selection criteria. 
The assumed energy resolution of a = 10 %/\/~E leads to the a m = 0.85 GeV resolution for 
the reconstructed invariant mass of the di-photon system originating from the H —>• 77 
decay. 
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Figure 4: Left: The multiplicity of hard-process isolated photons (left) and reconstructed 
isolated photons (right) for gg —>• H, H —>• 77 events with run = 125 GeV. 

2.5 Jets 

Clusters which have not been selected as associated with electrons or photons are smeared 
with Gaussian resolution parametrised in function Smearing: :forHadron (default: a = 
50 %/VE and 100 %/y/E). 
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Multiplicity 




Figure 5: Multiplicity of jets (top), A R cone distance between reconstructed jet and bary- 


center of particles (middle-left) and hard-process parton (middle right). The ratio of 
Pt t /pf^ irttcles (bottom-left) and p?* /Pt P parton (bottom-right) for gg —>• H,H —>• uu and 


mu = 125 GeV. 


8 




If the non-isolated muon falls into the cone of a cluster its 4-momenta is added to 
the cluster 4-momenta and the cluster direction is recalculated. The resulting clusters are 
classified as jets if their transverse momentum is greater than a given threshold (default: 
Pt > 15 GeV). They are removed from the OutputRecord.Clusters and stored in the 
common OutputRecord. Jets. 

The jets reconstruction efficiencies and di-jet mass resolution have been studied using 
control physics process of the gg —> H production with mu = 125 GeV and forcing the 
Higgs boson decay into specific partons, namely H -> uu, H -+ cc and H —»• bb. We 
used also process gg -+ H -+ rr for estimating tau-jet reconstruction efficiency (we forced 
tau-leptons decay into hadrons). For estimating reconstruction efficiency we consider 
only jets which have been reconstructed within the cone A R = 0.4 from the primary 
parton (particle) and we require that the primary parton (particle) has passed the same 
kinematical selection as required for reconstructed jets. 

2.5.1 Labeling 

Very important for the physics at LHC are jets originating from b-quarks (so called b-jets) 
which can be identified in the detector using b-tagging technique (vertex or soft-lepton 
tags). The package labels a jet as a b-jet if it is reconstructed within a limited rapidity 
range (default: \Ar/\ < 2.5) and if a b-quark of a transverse momenta (after FSR) above 
the threshold (default: pr > 5 GeV) is found within the cone (default: A R = 0.2) around 
the axis of reconstructed jet. The similar criteria are used for labeling the c-jets. 

Equivalently important are also jets originating from the hadronic r-decay (so called 
r-jets) which can be identified using dedicated algorithms. The package labels jet as 
a tau-jet if the hadronic decay product is relatively hard (default: pTffi had > 10 GeV), 
inside limited rapidity range (default: \r]\ < 2.5), dominates reconstructed jet transverse 
momenta (default: pffi had / > 0.9), and is within the cone (default: A Rj e t,r-had < 0.3) 
around the axis of a jet. 

Table 1 summarises the jets reconstruction+labeling efficiencies as obtained for the 
WH,H —* bb, cc, uu and gg —> H -+ rr events. Jets labeling is optional, can be switched 
off for b- and c-jets and/or separately for tau-jets (default: ON). 


Table 1: Efficiency for jet reconstruction+labeling for different types of initial partons with 
Pj? rton > 15 GeV (required p ^ 1 > 15 GeV). The rapidity coverage is limited to |? 7 | < 2.5. 
The A R cone = 0.4 is used for cluster reconstruction and AR cone = 0.2 is used for matching 
criteria. The gg —» H,H -+ bb, cc, wu and gg —> H —> rr processes were generated with 
mu = 125 GeV. In case of tau-jets only hadronic tau decays were generated. 


Parton type 

Reconstruction + Labeling 

u-quark 

95% 

b-quark 

81% 

c-quark 

87% 

tau-jet 

80% 
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Figure 6: Multiplicity of b-jets (top left) and hard-process b-quarks (top-right), A R 
cone distance between b-labeled jet and hard-process b-quark (bottom-left); the ratio of 
Pr t /PT Ph - quarks (bottom-right). Distributions are shown for gg —>• H,H —>• bb and 
run = 125 GeV. 
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Figure 7: Multiplicity of tau-jets (top left) and hadronic tau (top-right), A R cone distance 
between tau-labeled jet and hard-process tau-had (bottom-left); the ratio of p^ u ~ ]et f^ u ~ had 
(bottom-right). Distributions are shown for gg —>• H,H —>• rr and mjj = 125 GeV. 

2.5.2 Calibration 

The reconstructed jets four-momenta need to be corrected for the out-cone energy loss 
(cascade outside the jet-cone) and for the loss of the particles escaping detection (those 
below threshold at px = 0.5 GeV, neutrinos, invisible particles, muons outside acceptance 
range or below the observability threshold). Such correction, called calibration, can be 
performed on the statistical basis only. The single default calibration function (the same 
for any type of jets), as a function of transverse momenta of reconstructed jet, pjt 1 , is 
provided in the package. The calibration algorithm corrects jets four-momenta without 
altering their direction. Calibration can be switched off (default: ON). 

The quality of the calibration algorithm can be verified by monitoring the ratio of 
p T quark /p b T Jei , taking into account a hard-process quark which originates a given jet. 
One can note, (see Fig. 7, bottom plots) that the implemented calibration function has a 
tendency to undercalibrate b-jets while calibrates reasonably well u-jets. We can observe 
also rather large tail in the p^ 1 /p b T quark distribution, caused by the semileptonic b-quarks 
decays and larger spread of the cascading particles than in the case of light jets. 

As an effect of the hadronisation and cascading decays, the expected resolution for the 
resonance reconstruction in the hadronic channels will be much worse than in the leptonic 
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ones. The precision for reconstructing the peak position of the invariant mass of the di-jet 
system will relay on the precision of the calibration procedure. 

2.6 Missing transverse energy 

The missing transverse energy is calculated by summing up transverse momenta of iden¬ 
tified isolated photons, electrons and muons, of jets and clusters not accepted as jets and 
of non-isolated muons not added to any jet. Finally, the transverse energies deposited in 
cells not used for clusters reconstruction are also included in the total sum. Transverse 
energies deposited in unused cells are smeared with the same energy resolution function 
as for jets, and cells with deposited transverse energy below a given threshold (default: 0 
GeV) are excluded from the sum. From the calculation of the total sum E^ s the missing 
transverse energy is obtained, E™ iss = —Ej? s as well as the missing transverse momen¬ 
tum components p™ iss = —p° bs , p™ lss = — Py bs . The total calorimeter transverse energy, 
J2 Ej? l °, is calculated as the sum of all the above transverse energies except that of muons. 
Please note, that missing transverse energy is calculated from the energy balance before 
jets calibration is performed, as the possible out-cone energy loss is already taken into 
account by summing up energy deposition of unused cells/clusters. 

Fig. 8 shows the resolution of the transverse missing energy obtained for the di-jet 
events generated with the transverse momenta of the hard process above 17 GeV. 




Figure 8: The reconstructed missing transverse energy (left) and transverse momenta on 
neutrino inW^-eu events. 

2.7 Additional efficiencies 

The algorithms of the AcerDET package are not correcting for inefficiencies of photon, 
electron and muon reconstruction and identification. To be more realistic, one should 
apply a weighting factor of 70%-90% for each isolated lepton used in the analysis and of 
80% for each isolated photon. 

The package is also not correcting for tagging-efficiencies, namely the labeling proce¬ 
dure is not equivalent to the b- and tau-jet identification in the experiment. One can 
assume b-tagging efficiency of 60% per b-labeled jet with mistagging probability of 10% 
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for c-labeled jet and 1% for the light jet. For the tau-jets using efficiency of 50% (per 
tau-labelled jet) and 5 - 10 % mistagging probability for other jets could be a reasonable 
assumptions. One should be well aware that what proposed above represents quite crude 
estimates. 

The package is also not correcting for trigger efficiencies which has to be applied in 
addition if a given reconstructed object is foreseen to trigger an event. 

2.8 OUTPUT format 

By the end of the simulation and reconstruction algorithm, reconstructed entities: photons, 
electrons, muons, jets, transverse missing energy are available as collections of objects in 
the OutputRecord. Provided is algorithm to store this information in the ROOT [11] tree. 
However, the user might decide to use his preferred data-base technology for storing output 
from the simulation and reconstruction algorithms. 

3 Outlook 

We presented an update from Fortran to C++ of a package which can be useful for several 
phenomenological studies on the high px physics at LHC. As examples of applications 
some recent publications which have been using its fortran version AcerDET-1.0 [1] are 
papers eg. [14, 15, 16, 17] ( for the more complete list see link [2] ). 

It is not the aim of the package to represent in details performance neither of ATLAS 
nor of CMS detectors, nevertheless some global features of these would be reproduced 
well. We have presented some benchmark results in Section 2. In particular we believe 
that the analyses performed with the package for physics at LHC will be more realistic 
than parton-level studies alone. 

The package allows also for flexible adjusting of key parameters which characterise 
features of any detector be it for LHC or for future ILC or FFC experiments. 
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A General informations 


The simulation and reconstruction algorithm is executed by a call to the routine 
AcerDET. analyseRecord. The events which is going to be processed should be stored in 
the InputRecord. 

The convention for the particles status, mother-daughter relations, particles codes, etc. 
should be the same as in the HepMC standard. 

The input/output logical identifiers should be defined in conf igFileName. 

The package reads single input file acerdet .dat which contains parameters for events 
simulation and reconstruction and write out single control output file acerdet. out. If ini¬ 
tialised by the user in the main program, control histograms and ntuple with reconstructed 
events will be stored in the XXX . root file, where XXX stands for file which configures gen¬ 
erator used for generating events. In the distributed example the Pythia 8.1 is used. 

A.l Class AcerDET 

This is the main class, called to execute simulation and reconstruction algorithms. The 
following methods are implemented: 

• AcerDET contructor - initalisation, which should be called before the first event is 
processed. 

• AcerDET : : analyseRecord - simulation and reconstruction, should be called event 
by event. 

• AcerDET destructor - finalisation, should be called after last event is processed. 

TheAcerDET: : analyseRecord invokes other methods: analyse_Cell, analyse_Cluster, 
analyse_Muon, analyse_Electron, analyse_Photon, analyse_Jet, analyse_Mis, analyse_BJet, 
analyse_CJet, analyse_Tau, analyse_Calibration. The order in which methods are in¬ 
voked should not be changed. 

The AcerDET : : print Inf o() invokes priting information on the configurations used 
for individual methods and AcerDET : :printResults() summary information on the re- 
contructed objects. 

A.2 Interfaces to Event Generators 

The event, which is going to be processed, should be stored in the HepMC format. In the 
source code of the program provided are methods to traslate this format into internal 
event structure, vector of Particles stored in InputRecord and to interpret status and 
pdgld codes of the particles. 

The following internal types are introduced, based on the information available in the 
HepMC conversion of Pythia 8.1 event record. 

• ParticleType 

- PT_CJET - particle with pdgld code = -4, 4 

- PT_BJET - particle with pdgld code = -5, 5 

— PT_ELECTRON - particle with pdgld code = -11, 11 
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— PT_NEUTRINO_ELE - particle with pdgld code = -12, 12 

- PT_MUON - particle with pdgld code = -13, 13 

— PTJNTEUTRINOJMUON - particle with pdgld code = -14, 14 

- PT_TAU - particle with pdgld code = -15, 15 

— PT_NEUTRINO_TAU - particle with pdgld code = -16, 16 

- PT_PHOTON - particle with pdgld code = 22 

- PT_BOSON_Z - particle with pdgld code = 23 

— PT_BOSON_W - particle with pdgld code = -24, 24 

- PT_BOSON_H - particle with pdgld code = 25 

- PT_UNKNOWN - other 

• ParticleStatus 

— PSJ3EAM - corresponding HepMC: part— >is_beam() = true; 

— PS_FINAL - corresponding HepMC: part— >is_undecayed() = true; 

— PS_DECAYED - corresponding HepMC: part— >has_undecayed() = true; 

— PS_HISTORY - corresponding HepMC: part— >is_undecayed() = true; 

— PS_CASCADE_QUARK - corresponding HepMC: abs(part— >status()) ^,=30; 

- PS_HP_QUARK - corresponding HepMC: abs(part— >status()) = 20 - 29; 

- PS.NULL - other; 

While ParticleStatus and ParticleType become an attribute of the object Particle, 
the method AcerDet: : core: : isHardProcess is also provided to decide if a given particle 
is from hard process. This method (so far) relies on information of the particle origin, i.e. 
type of the mother particle, but more refined procedure can be impmeneted there as well. 

A.3 External calling sequence 

The following calling sequence should be provided in the main program (see demo.cpp): 

• initialise generator 

• initialise and configure AcerDET 

• initialise root tree and histogram manager 

• start even loop 

- generate new event 

- convert to HepMC format (if needed) 

- process with AcerDET 

- fill root tree with reconstructed objets 

• write tree and histograms into file, write information into log file. 
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Invoking method which fill root tree is optional. User however may decide to skip this 
part, or implement different output format. 

In general, if user do not wish, the dependencies on ROOT could be easily removed. 

It requires linking into executable the respective alternative library and providing very 
simple conversion package where the existing calls to ROOT functions for histograming are 
used to fill information into the user’s preferred ones. 

A.4 Execution sequence 

After the initialisation phase, the event by event execution sequence is a following one: 

• Map of the cells energy deposition is created: analyse_Cell. 

• Calorimetric clusters are reconstructed: analyse_Cluster. 

• Isolated and non-isolated muons are reconstructed: analyse_Muon. 

• Isolated electrons are reconstructed: analyse_Electron. 

• Isolated photons are reconstructed: analyse_Photon. 

• The remaining calorimetric clusters are identified as jets: analyse_Jet. 

• Missing transverse is calculated for the reconstructed event: analyse_Mis. 

• Optionally, if specified in the acerdet.dat file, algorithms for jets labeling and calibra¬ 
tion of jets are executed (default—ON): analyse_BJet, analyse_CJet, analyse_Tau, 
analyse_Calibration 

• Reconstructed event is stored in the final record: OutputRecord. 

The respective sequence of calls is executed in the AcerDET: : analyseRecord. 

A.5 Structure of the distribution version 

The distribution version consists of the source code of the AcerDET and example of the 
main program for execution with PYTHIA 8.1 generator. 
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Figure 9: Structure of AcerDET package. 

AcerDET is designed to be an universal tool. It does not depend on any external library 
(except C++ STL). In order to provide functionalities not built-in directly in AcerDET. it 
provides a set of interfaces. 

Instantiation of AcerDET requires providing instances of such interfaces. Because of 
separation between analyse logic and data structures AcerDET might work with a various 
number of external libraries. 

The AcerDET package provides a set of default implementations for its interfaces, 
grouped in AcerDET: :external subpackage containing: 

• Root_HistogramManager - a default implementation of IHistogramManager inter¬ 
face. Provides methods for storing weighted data in named histograms. For more 
details see AcerDET: : core : : IHistogramManager. 

• Pythia8_ParticleDataProviderFactory - a default implementation of 
IParticleDataProviderFactory interface. Represents a particle types database, 
providing particle charge, name, etc. For more details see 

AcerDET::core::ParticleDataProvider. 

• HepMC_InputConverter - static class converting HepMC event record to AcerDET in¬ 
ternal representation of input record. For more details see AcerDET : : io: : InputRecord. 

• RootJJTupleManager - wrapper for ROOT’S NTuple. Provides methods for storing 
and saving weighted events. The demo uses it like an external database. 

The acerdet. dat file (configuration file for AcerDET package) resides in main directory 
acerdet_dat. Source code resides in src which has one level deeper internal structure (see 
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README there). Interfaces to external libraries are in directory external The library 
will be created in src directory, provided is also example of the main program demo. cpp 
to execute AcerDET algorithms with PYTHIA 8.1 generator: 

• Struture of directory acerdet 


-rw-rw-r— 

i 

erichter 

erichter 

3703 

Jun 

10 

16:07 

acerdet.dat 

-rw-rw-r— 

i 

erichter 

erichter 

65 

Jun 

10 

16:07 

common.inc 

drwxrwxr-x 

2 

erichter 

erichter 

4096 

Jun 

12 

10:31 

conf 

-rw-rw-r— 

1 

erichter 

erichter 

2719 

Jun 

10 

17:25 

demo.cpp 

-rw-rw-r— 

1 

erichter 

erichter 

590 

Jun 

10 

16:07 

demoCreateConfig.cpp 

drwxrwxr-x 

3 

erichter 

erichter 

4096 

Jun 

12 

10:31 

external 

-rw-rw-r— 

1 

erichter 

erichter 

861 

Jun 

10 

16:08 

external.inc 

drwxrwxr-x 

8 

erichter 

erichter 

4096 

Jun 

12 

17:55 

.git 

-rw-rw-r— 

1 

erichter 

erichter 

74 

Jun 

10 

16:09 

.gitignore 

-rw-rw-r— 

1 

erichter 

erichter 

848 

Jun 

10 

16:07 

Makefile 

-rw-rw-r— 

1 

erichter 

erichter 

204 

Jun 

10 

16:08 

path.inc 

-rw-rw-r— 

1 

erichter 

erichter 

274 

Jun 

10 

16:07 

README 

-rw-rw-r— 

1 

erichter 

erichter 

69 

Jun 

10 

16:07 

setup.sh 

drwxrwxr-x 

6 

erichter 

erichter 

4096 

Jun 

12 

10:31 

src 


• Struture of directory acerdet/src (after compilation using ’make’) 


-rw-rw-r— 

i 

erichter 

erichter 

4078 

Jun 

10 

17:25 

AcerDET.cpp 

-rw-rw-r— 

i 

erichter 

erichter 

3994 

Jun 

10 

17:25 

AcerDET.h 

-rw-rw-r— 

i 

erichter 

erichter 

264412 

Jun 

10 

17:25 

AcerDET.o 

drwxrwxr-x 

2 

erichter 

erichter 

4096 

Jun 

12 

10:31 

analyse 

drwxrwxr-x 

2 

erichter 

erichter 

4096 

Jun 

10 

17:25 

conf 

drwxrwxr-x 

2 

erichter 

erichter 

4096 

Jun 

11 

19:14 

core 

-rw-rw-r— 

1 

erichter 

erichter 

99229 

Jun 

10 

16:08 

doxygen.config 

drwxrwxr-x 

2 

erichter 

erichter 

4096 

Jun 

10 

17:25 

io 

-rw-rw-r— 

1 

erichter 

erichter 

5538158 

Jun 

12 

10:31 

libAcerDET.a 

-rwxrwxr-x 

1 

erichter 

erichter 

3803523 

Jun 

12 

10:31 

libAcerDET.so 

-rw-rw-r— 

1 

erichter 

erichter 

1689 

Jun 

10 

16:07 

Makefile 

-rw-rw-r— 

1 

erichter 

erichter 

132 

Jun 

10 

16:08 

path.inc 

-rw-rw-r— 

1 

erichter 

erichter 

252 

Jun 

10 

16:07 

README 


To execute the package following actions should be taken: 

• In subdirectory acerdet edit files external. inc and path.inc. Then type: make . 
It will compile the code and create the src/libAcerDET. a and src/libAcerDET. so. 
It will create also demo executable demo. exe 

• To execute, type 

./demo.exe conf/XXX.conf nEvent 

where conf/XXX. conf stands for file configuring generated process and nEvent is a 
integer number of required events to be generated. 

• Output will be stored in hie conf/XXX. conf .root 


A.6 Development 

Regarding modifications in AcerDET algorithms, package users are able to develop own 
extensions to our main algorithm. As it was shown in previous section AcerDET provides 
a set of interfaces to implement in order to use the package. AcerDET: : external sub¬ 
package contains default implementations of those interfaces, but user may develop his 
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own implementations. For instance if someone wants to use other histogram manager, it 
is enough to create new class inheriting from AcerDET: : core: : IHistogramManager and 
pass it’s instance when create new AcerDET instance (as it is shown in demo.cpp). 


const std::string configFileName = "acerdet.dat"; 

Configuration configuration = Configuration::fromFile( configFileName ); 

IParticleDataProviderFactory *pdpFactory = 

new external::Pythia8_ParticleDataProviderFactory(); 

IHistogramManager *histoManager = 

new external::Root_HistogramManager(); // HERE USE OWN INSTANCE 
external::Root_NTupleManager nTuple; 

AcerDET acerDet( 
configuration, 
pdpFactory, 
histoManager 
); 


Analoguously user may use own instance of IParticleDataProviderFactory. ROOT 
NTuple is used in demo like na external database, so using other databases is possible in 
similiar way. Distributed demo version uses HepMC as input data standard, while AcerDET 
uses its own input format. Default conversion algoritm is placed in 

AcerDET: :external: :HepMC_InputConverter. In order to use other input format, simil¬ 
iar converter to AcerDET input format is required. 
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B List of input parameters: acerdet.dat file 


Flag.HistogramlD 

10000 

....id for histograms 




Flag.Smearing 

1 

....smearing on=l, off=0 




Flag.B-Field 

1 

. . ..B-field on=l, off=0 




Flag.SUSY_LSP_Particle 

66 

....code for SUSY LSP particle 




Flag.BC-JetsLabeling 

1 

....b- and c-jets labeling on=l, off=0 




Flag.Tau-JetsLabeling 

1 

....tau-jets labeling on=l, off=0 




Flag.JetCalibration 

1 

....jet calibration on=l, off=0 




Cell.RapidityCoverage 

5.000 

....rapidity coverage 




Cell.MinpT 

0.500 

....min p_T for B-field 




Cell.MinEt 

0.000 

....min E_T for cell 




Cell.EtaTransition 

3.200 

....eta transition in cells granularity 




Cell.GranularityEta 

0.100 

....granularity in eta (within Cell.EtaTransition), 

2x 

outside 

Cell.GranularityPhi 

0.100 

....granularity in phi (within Cell.EtaTransition), 

2x 

outside 

Cluster.MinEt 

5.000 

....minimum E_T for cluster 




Cluster.ConeR 

0.400 

....cone R for clustering 




Cluster.RapidityCoverage 

5.000 

....rapidity coverage 




Cluster.MinEtlnit 

1.500 

....min E_T for cluster initiator 




Muon.MinMomenta 

6.000 

....minimum muon-momenta to be detected 




Muon.MaxEta 

2.500 

....maximum muon eta to be detected 




Muon.MinlsolRlj 

0.400 

....min R_lj for muon-isolation 




Muon.ConeR 

0.200 

....R_cone for energy deposition 




Muon.MaxEnergy 

10.000 

....max energy deposition for isol 




Photon.MinMomenta 

5.000 

....minimum photon-momenta to be isol 




Photon.MaxEta 

2.500 

....maximum photon eta to be isol 




Photon.MinJetsRlj 

0.150 

....min R_lj for photon-jet 




Photon.MinlsolRlj 

0.400 

....min R_lj for photon-isolation 




Photon.ConeR 

0.200 

....R_cone for energy deposition 




Photon.MaxEnergy 

10.000 

....max energy deposition for isol 




Electron.MinMomenta 

5.000 

....minimum electron-momenta to be isol 




Electron.MaxEta 

2.500 

....maximum electron eta to be isol 




Electron.MinJetsRlj 

0.150 

....min R_lj for electron-jet 




Electron.MinlsolRlj 

0.400 

....min R_lj for electron-isolation 




Electron.ConeR 

0.200 

....R_cone for energy deposition 




Electron.MaxEnergy 

10.000 

....max energy deposition for isol 




Jets.MinEnergy 

10.000 

....jets energy_min threshold 




Jets.RapidityCoverage 

5.000 

....rapidity coverage for jets 




BJets.MinMomenta 

5.000 

....minimum b-quark pT (after FSR) momenta 

for b-jet 

label 

BJets.MaxEta 

2.500 

....maximum b-quark eta for b-jet label 




BJets.MaxRbj 

0.200 

....max R_bj for b-jet label 




CJets.MinMomenta 

5.000 

....minimum c-quark pT (after FSR) momenta 

for c-jet 

label 

CJets.MaxEta 

2.500 

....maximum c-quark eta for c-jet label 




CJets.MaxRcj 

0.200 

....max R_cj for c-jet label 




Tau.MinpT 

10.000 

....minimum tau-had pT for tau-jet label 




Tau.MaxEta 

2.500 

....maximum tau-eta for tau-jet label 




Tau.MinR 

0.300 

....max R_tauj for tau-jet 




Tau.MaxR 

0.900 

....max R_tauj for tau-jet 




Misc.MinEt 

0.000 

....min E_T for energy in cell to count unused cell 
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C Reconstructed entities 


Reconstructed entities are stored as vectors of ObjectData (for isolated photon, muons, 
electrons), as CellData (for cells), as JetData (for jets) and MisData for missing energy. 
Objects in each category of ObjectData and JetData are ordered with the decreasing 
transverse momenta. 

C.l ObjectData 

For ObjectData stored is information of the pdg_id, so the information of the charge 
is preserved, eta, phi and pT which allows to recalculate full kinematics assuming that 
mass is equal to zero, and the internal flag alreadyUsed needed for calculating missing 
tranverse energy for the event. 

C.2 JetData 

For JetData stored is information on the type which flags if a given jet is b-quark tagged, 
c-quark tagged or matched with hadronic-tau decays. The angular kinematic variables 

eta, phi and eta_rec, phi_rec and pT. 

C.3 MisData 

For MisData stored is information on reconstructed momenta both transverse compo¬ 
nents of electrons, muons, photons, clusters and jets, PXSUM, PYSUM, on reconstructed 
total tranverse components i.e. after adding non-clusterd cells PXREC, PYREC, sum of 
transverse momenta components, the reconstructed transverse momenta SUMET, the trans¬ 
verse components of momenta of neutrinos and invisible particles (as defined by the user) 
in the event PXNUES, PYNUES and finally sum of the total transverse energy components 
deposited in the calorimeter PXXCALO, PYYCALO. 

D Content of the ACDTree 

Reconstructed objects and, for convenience of the analysis, also partial information on 
generated event are stored in the format of ROOT tree, the ACDTree. 

• Generator info: Event weight and identifier for generated process. 

• Generated particles: stored are those participating in the hard-scattering, which is 
defined using isHardProcess method and flags: PS_Beam, PS_HISTORY. Stored are 
4-momenta, and pdg_id code. 

• Isolated photons, electrons, muons: stored are 4-momenta, and pdg_id code. 

• Jets: stored are 4-momenta, the pdg_id code indicate if jet was labeled as -batteg, 
c-tagged or hadronic-decay of tau lepton. 

• Missing energy: stored are transverse momenta components of missing energy, in¬ 
visible energy and energy reconstructed in the calorimeter. 
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Additional information, if required, can be also extracted from the internal variables of 
simulation and reconstruction algorithms: eg. on used/unused cells, non-isolated muons, 
unused clusters, etc. and added to the ACDTree. 


E Parametrisation for energy/momenta resolution 

E.l Smearing::forPhoton 

The parametrisation for photon energy resolution assumes only energy dependence and 
the Gaussian smearing with 10 %/^/E^ resolution. 

E° meared = E^ ue • (1 + r n • °' 10 ). (1) 

yjEtfue 

Where r n is the random number generated according to the Gaussian distribution (normal 
distribution from Standard C++ library). All 4-monrenta components of the photon are 
smeared with the same resolution so the direction of the photon is not altered. 


E.2 Smearing:: forElectron 

The parametrisation for electron energy resolution assumes only energy dependence and 
the Gaussian smearing with 12%/y r E^ resolution. 


■^smeared 


= E. 


true 


■ (1 + r n 


0.12 

yjEl rUe ' 


( 2 ) 


Where r n is the random number generated according to the Gaussian distribution (normal 
distribution from Standard C++ library). All 4-monrenta components of the electron are 
smeared with the same resolution so the direction of the electron is not altered. 


E.3 Smearing: :forHadron 


The parametrisation for clusters momenta resolution assumes only energy dependence and 
Gaussian smearing with 50%/V E c i u or 100 %/y/Edu resolution. The transition region, 
| rf :lu | = CALOTH is the same as the transition of the cells granularity. 

For \rj du \ < CALOTH: 


rpsmeared _ jptrue 

&clu — &clu 


■ (1 + r n 


0.50 . 

/ zptrue 
V ^clu 


( 3 ) 


For \r] clu \ > CALOTH: 


Tpsmeared _ rptrue 

i^du ~ &clu 


(1 + r n 


1.00 ^ 

/ jT'true 
V ^clu 


( 4 ) 


Where r n is the random number generated according to the Gaussian distribution (normal 
distribution from Standard C++ library). All 4-monrenta components of the muon are 
smeared with the same resolution so the direction of the muon is not altered. 


24 



E.4 Smearing::forMuon 


The parametrisation for muon momenta resolution assumes only transverse momenta de¬ 
pendence and Gaussian smearing with 0.0005 • pT resolution. 

pT smeared = pT j[ ue /( 1 + r n • 0.0005 • pT^ ue ). (5) 

Where r n is the random number generated according to the Gaussian distribution (normal 
distribution from Standard C++ library). All 4-momenta components of the muon are 
smeared with the same resolution so the direction of the muon is not altered. 

E.5 Particle::foldPhi 

The effect of the magnetic field is included by simple shifting the cf> position of the charged 
particles respectively to its transverse momenta, parametrised as following: 

S<j> = 0.5 /pT part . (6) 

The 5(f) is calculated in radians and pT part is given in GeV. The sign of the 5<f> is the same 
as the sign of the particle charge. Charged particles with pT part < 0.5 GeV are assumed 
to be looping in the detector and not depositing energy in the calorimeter. 

F Some formulas 

We work with the assumptions of massless reconstructed objets. The following relations 
were used to translate between (pT,ip(j)) coordinates and four-momenta (p x ,p y ,p Zl E). 


Px = PT ■ COs((f> ) 

(7) 

Py = PT ■ sin{4>) 

(8) 

Pz= PT- cosh(rj) 

(9) 

E = pt ■ sinh(rj) 

(10) 

PT = \jpl+P 2 y 

(11) 

„ Jpt + P 2 z + \Pz\ 
sign( In ,p z ) 

PT 

(12) 

(f> = asinh{py/^p 2 x +p 2 y ) 

(13) 
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G Output content 

G.l Output ACDTree content: acerdet.ntup file 

The content of the ntuple is exactly as described in the previous section. As an auxiliary 
info, within the stored information, could be provided code on the generated process, 

IDPROC. 


3|e3)e^e3|e^e3)o|e3|e^e3f;^e3f;^e3f;9|e3f;3|e3)o|e3f;^e3f:3f;3|e3|e^e^e3)o|c^c3|e^c3|e3|e^e3|c3|c^c3|e^c^c3|c^:^c3|e^c3|e3|c^c3|c3|e^e3| c: t :: t : H( : 4 c: 4 c H c: t : H‘3k^ e H e3 l c: 4 c H‘ : t :: t : H c: 4 c: t : H‘ : 4 e H c H( 3 l c H‘ 

*Tree :ACDTree : ACDTree * 

*Entries : 100000 : Total = 176704219 bytes File Size = 83438909 * 

* : : Tree compression factor =2.12 * 

3|e^e^e^e3|e3f;3|e^e^e3f;3|e3|e^e3f;^e^e^e^e3|e^e3|e^e^e3|e^e3|e3|e3)o|e^c3|e^c3|e3|e^c3|e3|c^c3|e^c^c3|c^c^:3|c^c3|e3|c^:3|c3|e^e3|c^:^c3|e^:^:3|c^:3|e3|e^e3|o|c^c3|c^:^c3|c^c^c3|e^e3|c3|c^c3|e 


*Br 0 

*Entries 

*Baskets 

ProcessID : ProcessID/I * 

100000 : Total Size= 785597 bytes File Size = 401027 * 

13 : Basket Size= 32000 bytes Compression= 1.96 * 

*Br 1 

*Entries 

*Baskets 

part_n : part_n/I * 

100000 : Total Size= 785546 bytes File Size = 400988 * 

13 : Basket Size= 32000 bytes Compression= 1.96 * 

*Br 2 

*Entries 

*Baskets 

part_pdgld : vector<int> * 

100000 : Total Size= 14579602 bytes File Size = 7086584 * 

235 : Basket Size= 32000 bytes Compression= 2.06 * 

*Br 3 

*Entries 

*Baskets 

part_mothej 

100000 

235 

-_pdgld : vector<int> * 

Total Size= 14581275 bytes File Size = 7088229 * 

Basket Size= 32000 bytes Compression= 2.06 * 

*Br 4 

*Entries 

*Baskets 

part_px 

100000 

235 

vector<float> * 
Total Size= 14578885 bytes File Size = 7085879 * 
Basket Size= 32000 bytes Compression= 2.06 * 

*Br 5 

*Entries 

*Baskets 

part_py 

100000 

235 

vector<float> * 
Total Size= 14578885 bytes File Size = 7085879 * 
Basket Size= 32000 bytes Compression= 2.06 * 

*Br 6 

*Entries 

*Baskets 

part_pz 

100000 

235 

vector<float> * 
Total Size= 14578885 bytes File Size = 7085879 * 
Basket Size= 32000 bytes Compression= 2.06 * 

*Br 7 

*Entries 

*Baskets 

part_E 

100000 

235 

vector<float> * 
Total Size= 14578646 bytes File Size = 7085644 * 
Basket Size= 32000 bytes Compression= 2.06 * 

*Br 8 

*Entries 

*Baskets 

ele_n 

100000 

13 

ele_n/I * 
Total Size= 785529 bytes File Size = 400975 * 
Basket Size= 32000 bytes Compression= 1.96 * 

*Br 9 

^Entries 

*Baskets 

ele_pdgld 

100000 

57 

vector<int> * 
Total Size= 3204969 bytes File Size = 1411515 * 
Basket Size= 32000 bytes Compression= 2.27 * 

*Br 10 

*Entries 

*Baskets 

ele_px 

100000 

57 

vector<float> * 
Total Size= 3204786 bytes File Size = 1411344 * 
Basket Size= 32000 bytes Compression= 2.27 * 

*Br 11 

*Entries 

*Baskets 

ele_py 

100000 

57 

vector<float> * 
Total Size= 3204786 bytes File Size = 1411344 * 
Basket Size= 32000 bytes Compression= 2.27 * 

*Br 12 

^Entries 

*Baskets 

ele_pz 

100000 

57 

vector<float> * 
Total Size= 3204786 bytes File Size = 1411344 * 
Basket Size= 32000 bytes Compression= 2.27 * 
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*Br 13 

*Entries 

*Baskets 

ele_E 

100000 

57 

vector<float> * 
Total Size= 3204725 bytes File Size = 1411287 * 
Basket Size= 32000 bytes Compression 2.27 * 

*Br 14 

*Entries 

*Baskets 

muo_n 

100000 

13 

muo_n/I * 
Total Size= 785529 bytes File Size = 400975 * 
Basket Size= 32000 bytes Compression 1.96 * 

*Br 15 

^Entries 

*Baskets 

muo_pdgId 

100000 

57 

vector<int> * 
Total Size= 3204997 bytes File Size = 1411543 * 
Basket Size= 32000 bytes Compression 2.27 * 

*Br 16 

^Entries 

*Baskets 

muo_px 

100000 

57 

vector<float> * 
Total Size= 3204814 bytes File Size = 1411372 * 
Basket Size= 32000 bytes Compression 2.27 * 

*Br 17 

^Entries 

*Baskets 

muo_py 

100000 

57 

vector<float> * 
Total Size= 3204814 bytes File Size = 1411372 * 
Basket Size= 32000 bytes Compression 2.27 * 

*Br 18 

^Entries 

*Baskets 

muo_pz 

100000 

57 

vector<float> * 
Total Size= 3204814 bytes File Size = 1411372 * 
Basket Size= 32000 bytes Compression 2.27 * 

*Br 19 

*Entries 

*Baskets 

muo_E 

100000 

57 

vector<float> * 
Total Size= 3204753 bytes File Size = 1411315 * 
Basket Size= 32000 bytes Compression 2.27 * 

*Br 20 

*Entries 

*Baskets 

pho_n 

100000 

13 

pho_n/I * 
Total Size= 785529 bytes File Size = 400975 * 
Basket Size= 32000 bytes Compression 1.96 * 

*Br 21 

*Entries 

*Baskets 

pho_pdgId 

100000 

78 

vector<int> * 
Total Size= 4558940 bytes File Size = 2093066 * 
Basket Size= 32000 bytes Compression 2.18 * 

*Br 22 

*Entries 

*Baskets 

pho_px 

100000 

78 

vector<float> * 
Total Size= 4558694 bytes File Size = 2092832 * 
Basket Size= 32000 bytes Compression 2.18 * 

*Br 23 

*Entries 

*Baskets 

pho_py 

100000 

78 

vector<float> * 
Total Size= 4558694 bytes File Size = 2092832 * 
Basket Size= 32000 bytes Compression 2.18 * 

*Br 24 

^Entries 

*Baskets 

pho_pz 

100000 

78 

vector<float> * 
Total Size= 4558694 bytes File Size = 2092832 * 
Basket Size= 32000 bytes Compression 2.18 * 

*Br 25 

^Entries 

*Baskets 

pho_E 

100000 

78 

vector<float> * 
Total Size= 4558612 bytes File Size = 2092754 * 
Basket Size= 32000 bytes Compression 2.18 * 

*Br 26 

*Entries 

*Baskets 

jet _n 

100000 

13 

jet_n/I * 
Total Size= 785529 bytes File Size = 400975 * 
Basket Size= 32000 bytes Compression 1.96 * 

*Br 27 

*Entries 

*Baskets 

iet_pdgld 

100000 

85 

vector<int> * 
Total Size= 4991805 bytes File Size = 2301791 * 
Basket Size= 32000 bytes Compression 2.17 * 

*Br 28 

*Entries 

*Baskets 

jet_px 

100000 

85 

vector<float> * 
Total Size= 4991538 bytes File Size = 2301536 * 
Basket Size= 32000 bytes Compression 2.17 * 
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*Br 29 

*Entries 

*Baskets 

jet-py 

100000 

85 

vector<float> * 
Total Size= 4991538 bytes File Size = 2301536 * 
Basket Size= 32000 bytes Compression 2.17 * 

*Br 30 

^Entries 

*Baskets 

jet_pz 

100000 

85 

vector<float> * 
Total Size= 4991538 bytes File Size = 2301536 * 
Basket Size= 32000 bytes Compression 2.17 * 

*Br 31 

*Entries 

*Baskets 

jet_E 

100000 

85 

vector<float> * 
Total Size= 4991449 bytes File Size = 2301451 * 
Basket Size= 32000 bytes Compression 2.17 * 

*Br 32 

^Entries 

*Baskets 

pxmiss 

100000 

13 

pxmiss/F * 
Total Size= 785546 bytes File Size = 400988 * 
Basket Size= 32000 bytes Compression 1.96 * 

*Br 33 

^Entries 

*Baskets 

pymiss 

100000 

13 

pymiss/F * 
Total Size= 785546 bytes File Size = 400988 * 
Basket Size= 32000 bytes Compression 1.96 * 

*Br 34 

*Entries 

*Baskets 

pxnue 

100000 

13 

pxnue/F * 
Total Size= 785529 bytes File Size = 400975 * 
Basket Size= 32000 bytes Compression= 1.96 * 

*Br 35 

^Entries 

*Baskets 

pynue 

100000 

13 

pynue/F * 
Total Size= 785529 bytes File Size = 400975 * 
Basket Size= 32000 bytes Compression= 1.96 * 

*Br 36 

^Entries 

*Baskets 

pxcalo 

100000 

13 

pxcalo/F * 
Total Size= 785546 bytes File Size = 400988 * 
Basket Size= 32000 bytes Compression= 1.96 * 

*Br 37 

*Entries 

*Baskets 

pycalo 

100000 

13 

pycalo/F * 
Total Size= 785546 bytes File Size = 400988 * 
Basket Size= 32000 bytes Compression= 1.96 * 


* . * 

r 
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G.2 Control printout 

^0|c^c3|e^c3|e^c3|e^c3|c^c3|e^c3|c^c3|e^c3|e^c3|e^c3|c3|c^c3|c^c^c3|c^c^c3|e^c3|c3|e^c3|e3|e^c3|c^c^c3|c^c^c3|c^c3|e3|c^c3|c3|c^c3|e^c^c3|c^c^c 


* AcerDET, version: 2.0 * 

* Released at: 30.06.2015 * 

* * 

* Simplied event simulation and reconstruction package * 

* * 

* by E. Richter-Was (Institute of Physics) * 

* and P. Mikos (Theoretical Computer Science) * 

* Jagiellonian University, Cracow, Poland * 


9|0|c^c3|e^c3|e^c3|e^c3|c3|c3|e^c3|c^c3|c^c3|c^c3|e3|c3|c3|c^c3|c^c^c3|c3|c^c3|e^c3|e3|e^c3|0|e^c3|c^c^c3|e^c^c3|c3|c3|e3|e^c3|e3|c^c3|e^c^c3|c^c^c 

Initial configuration: 

3fe^e^e3|e3|e3f;^e3f;^e^;^e3f;^e3f;^e3f;^e^;^e3f;3|e4;3|c^c3|e^c^c3|c^c^c3|e^c3|c3|e^c3|e3|e 


* * 

a|e s|e s|e s|e s)c 3|e s|e s|e j)c s|e j)c s|e s)e s)e s|e s|t s)e 

* *** analyse::Cell *** * 

a|c :|e :|e :|e :|e :|e :|e :|e 3|e :|e :|e :|e :|e sje 

* * 


^(3f;^e^e^e3f;3|e^e3|e3)o|e3f;3|e3f;^e^e3|e3f;3|e^e^e4;3|e^:3|c^c^:3|c^c^c3|c^c3|e3|e^c3|c3|c 

clusters definition ... 
eta coverage 5.000000 
E_T_min cell thresh 0.000000 
eta gran, transition 3.200000 
gran in eta(central) 0.100000 
gran in phi(central) 0.100000 

B field apply .... 

B-field on 

p_T min non looping 0.500000 

^(3f;^e^e^e^e^e3f;^e^;^e3f:3|e3f;^e3f:3|e^c^e3f;^e^o|c3|c3|e^c^c3|c3|c^c3|e^c3|e3|e^c3|e3|e^c3|c^c 


* * 

* *** analyse::Cluster *** * 

* * 


3fe:)e^e3|e^e3f;9|e3|e^e3f;^e3f;^e3f;9|e3f;^e3f;9|e3f;^e3f;3|e^c3|c^c^c3|e^e^e3|c^c3|e3|e^e3|c3|c^:3|e^c 

clusters definition ... 

E_T_min cluster 5.000000 

E_T_min cell initia 1.500000 

R cone 0.400000 

eta coverage 5.000000 

eta gran, transition 3.200000 

^C^e3|e^e^c3|e3|e3f;3|e^e^e3f;^e3f;^e3f;^e3|e^e^e^e^0|c3|c3|e^c^c3|c3|c^c3|e^c3|c3|e^c3|e3|e 


* * 

* *** analyse::Muon *** * 

3)c 3|c3|e3|c3|e3(e3|e3|e3|e3(c3|e3|c3|e3(c3|e3|e3(e3|c3|e3|e3|e3|e3|e3|e3|c3|e 

* * 


3fc3tc:(e3(c3fc3te3fc3tc3fe3tc3fe3(c3fc3tc3fc3tc3fe3fc3fe3(c3fc3fe3|e3|e3|e3(c3|c3|e3|e3|e3|e3|e3|e3|e3|e3|e3|e 

... muon isolation ... 

min. muon p_T 6.000000 

max. muon eta 2.500000 

min R_lj for isolation 0.400000 

R for energy deposit 0.200000 

max E_dep for isolation 10.000000 

smearing on 
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3|o|e^c3|e^c3|e^c3|e^c3|c^c3|e^c3|c^c3|e^c3|c^c3|e^c3|e3|e^c3|c^c^c3|e^c^c3|e^c3|c3|c^c3|e3|e^c3|c^c^c 

* * 

* *** analyse::Electron *** * 

3)c ^c s|e j)c :|e s(e 3|e ^c +; ^c + ^e s|c s|e j)c 3|e s|e j)c s|e s)c s|e j)c s|e s|c s)e 3|e s|e 

* * 
3|o|e^e3|c^c3|e^c3|e^c3|e^c3|e^c3|e^c3|e^c3|e^c3|c^c3|e3|e^c3|c^:^:3|e^c^c3|e^c3|e3|e^c3|o|c^c3|e^c^: 

... electron isolation ... 
min. lepton p_T 5.000000 
max. lepton eta 2.500000 
max R_ej for ele-clu 0.150000 
min R_lj for isolation 0.400000 
R for energy deposit 0.200000 
max E_dep for isolation 10.000000 
smearing on 

3|e3|c^c3|e^c3|e^c3|c^:3|c^:3|c^c3|e^c3|c^:3|c^c3|c^c3|e3|e^c3|c^:^:3|c^c^c3|e^c3|e3|e^c3|e3|e^c3|e 

* * 

s|e s|e j)c :|e s(e +: s|c s|e s|e s|e j)c s|e j)c s|e s)c s|t 3|e s|e s|c j)c s|e 3|e 

* *** analyse::Photon *** * 

% :+: j(c :+: :+: 3|e 

* * 
3fe3|e3|c3|e3|e3|e3|c3|e3(c3|e3(c3|e3|c3|e3(e3|e3(c3|e3|e3|e3|e3|e3|e3(:3|e3(c3|e3|c3|:3|e3|e3|c3|e3|e3|e3|e3|e3|c3|e 

... photon isolation ... 
min. photon p_T 5.000000 
max. photon eta 2.500000 
max R_gam-clust 0.150000 
min R_isol 0.400000 
R for energy deposit 0.200000 
max E_dep for isolation 10.000000 
smearing on 

3|o|e^c3|e^c3|e^c3|e^c3|e^c3|e^c3|e^c3|c^c3|c^c3|c^c3|e3|e^c3|c^:^:3|e^e^e3|e^c3|e3|e^e3|e 

* * 

s|c :|e :|e :|e 3|e :|c :|e 3|c 3|e :|e :|e :|e 3|c 

* *** analyse::Jet *** * 

:+c 3|e 

* * 
3|e3|c^c3|e^c3|c^c3|e^c3|c^c3|e^c3|c^c3|c^c3|c^c3|e^c3|e3|c^c3|c^c^:3|c^c^c3|e^c3|c3|c^c3|c 

clusters definition .... 

R cone 0.400000 
jets definition .... 

E_T_jets [GeV] 10.000000 
eta coverage jets 5.000000 
smearing on 
B-field on 
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3|o|e^c3|e^c3|e^c3|e^c3|c^c3|e^c3|c^c3|e^c3|c^c3|e^c3|e3|e^c3|c^c^c3|e^c^c3|e^c3|e3|e^c3|e 


* * 

:+c 3|e 

* *** analyse::Mis *** * 

3)c ^c :|e ^c :|e ^c :|e ^c :|e ^c :|c ^e :|c ^c :|e ^c 3|e :|c 3|e ^c :|e ^c sje 

* * 


3|o|e^e3|c^c3|e^c3|c^c3|e^c3|e^c3|e^c3|c^c3|c^c3|c^:3|e3|e^c3|c^:^:3|e^e^e3|e^c3|e3|e^e3|e 

muon coverage 
min. muon p_T 6.000000 
max. muon eta 2.500000 
unused cells . . . 
smearing on 

cells threshold 10.000000 
invisible particles ... 

KF code for invis 66 

3fc3|e3|c3|e3|c3|e3|c3|e3(c3|e3(c3|e3|c3|e3(c3|e3(e3|e3(c3|e3(c3|e3|e3|e3|e3(c3(c3|e3(e3|c3|e3(c3|e3|e3|c3|e3|e 


* * 

s|e :|e :|c :|e :|e :|e :|c :|e :|e :|e :|c 3|e :|e :|e s|c 

* *** analyse::BJet *** * 

% s|e s|e :|e s|e s)c + s|e s|e j)c s|e j)c j|e j|e j)c 3|e s|t s|e 

* * 


^(3|c^c3|e^c3|e^c3|c^c3|c^:3|c^c3|e^:3|c^:3|c^c3|c^c3|e3|e^:3|c^:^:3|c^c^e3|c^c3|e3|e^e3|c3|c 

... jets labeling ... 
labeling on/off on 
bjets ... 

min b-quark p_T 5.000000 
max b-quark eta 2.500000 
max R_bj for b-jets 0.200000 

3|e3|c^c3|e^c3|c^c3|e^c3|c^c3|e^c3|c^c3|e^c3|c^c3|e^c3|c3|c^c3|c^c^c3|c^c^c3|e^c3|c3|e^c3|e3|e 


* * 

sf: s|e s|e s|e s|e s|e + s|e j)c s|e j)c s|e s|e s)e 3|e s|t s|e 

* *** analyse::CJet *** * 

3|e :|e :|e :|c :|e :|e :|e :|e 3|e :|e :|e :|e :|e 3|c 

* * 


^(3|e^c3|c^c3|e^c3|c^c3|c^:3|c^:3|e^:3|c^:3|c^c3|c^c3|e3|e^:3|c^:^:3|c^c^c3|e^c3|e3|e^e3|c3|e 

... jets labeling ... 
labeling on/off on 
cjets ... 

min c-quark p_T 5.000000 
max c-quark eta 2.500000 
max R_cj for c-jets 0.200000 

^(3|c^c3|e^c3|c^c3|e^c3|c^c3|c^c3|c^c3|e^c3|c^c3|e^c3|c3|c^c3|c^c^c3|c^c^c3|e^c3|c3|e^c3|e 


* * 

s|c :|e :|e :|e :|e :|e :|e :|e :|e :|e :|e :|e 3|e 

* *** analyse::Tau *** * 

s|e s(c s|e s)c :|e s)c s|e s)c +; s|e s|e +: s)c s|t j)c j|c j)e s|e j)c j|e + s|c 3|c 

* * 


3|0|e^c3|e^c3|c3|c3|e^c3|c^c3|e^c3|c^c3|c^c3|c^c3|e^c3|c3|c^c3|c^c^:3|c^c^c3|e^c3|c3|c^c3|c 

... tau-jets labeling ... 
labeling on/off on 
tau-jets ... 

min tau-had p_T 10.000000 
max tau-had eta 2.500000 
max R_tauj for tau-jets 0.300000 
tau-had frac. of jet 0.900000 



9|o|e^c3|e^c3|e^c3|e^c3|c^c3|e^c3|c^c3|e^c3|c^c3|e^c3|e3|e^c3|c3|c^c3|e^c^c3|c^c3|e3|e^c3|e3|e^c3|c^c^c3|c^c^c 


* * 
* *** analyse::Calibration *** * 

s|e ^c :|e ^c :|e ^c :|c ^c :|e :|c ^c :|e ^c :|c :|e :|e 3|e ^c :|e ^c :|e ^c :|c :|c 3|e s|e ^c 


jets calibration .... 
calibration on 

G.3 Control histograms: XXX.conf.root file 

Below is the list of control histograms which monitors performance of the simulation/reconstruction. 
The are store together with ACDTree tuple in file XXX.conf.root. 

OBJ: TTree ACDTree ACDTree : 0 at: 0x8190c80 

KEY: TH1F hist10000;1 Cell: multiplicity 

KEY: TH1F histl0101;l Cluster: multiplicity 

KEY: TH1F histlOllljl Cluster: delta eta cluster barycentre 

KEY: TH1F hist10112;1 Cluster: delta phi cluster barycentre 

KEY: TH1F hist10113;1 Cluster: delta r cluster barycentre 

KEY: TH1F hist10123;1 Cluster: delta r cluster HPparton 

KEY: TH1F hist10114;1 Cluster: pTclu/SumpT particle 

KEY: TH1F hist10124;1 Cluster: pTclu/SumpT HP parton 

KEY: TH1F hist10210;1 Muon: muon multiplicity N0IS0LATED 

KEY: TH1F histl0211;l Muon: muon multiplicity ISOLATED 

KEY: TH1F histl0221;l Muon: muon multiplicity HP 

KEY: TH1F histl0231;l Muon: muon multiplicity HP+IS0L 

KEY: TH1F histl0311;l Electron: multiplicity ISOLATED 

KEY: TH1F histl0321;l Electron: multiplicity HP 

KEY: TH1F histl0331;l Electron: multiplicity HP+IS0L 

KEY: TH1F histl0411;l Photon: photon multiplicity ISOLATED 

KEY: TH1F histl0421;l Photon: photon multiplicity HP 

KEY: TH1F histl0431;l Photon: photon multiplicity HP+IS0L 

KEY: TH1F histl0501;l Jet: multiplicity 

KEY: TH1F histl0511;l Jet: delta phi jet-barycentre 

KEY: TH1F histl0512;l Jet: delta eta jet-barycentre 

KEY: TH1F histl0513;l Jet: delta r jet-barycentre 

KEY: TH1F histl0514;l Jet: delta r jet- HP parton 

KEY: TH1F histl0523;l Jet: pTjet/pT particles in Rcone 

KEY: TH1F hist10524;1 Jet: pTjet/pT HP parton 

KEY: TH1F histl0611;l Mis: reconstructed p_T 

KEY: TH1F hist10612;1 Mis: reconstructed p_T + cells 

KEY: TH1F hist10613;1 Mis: reconstructed pTmiss 

KEY: TH1F histl0621;l Mis: true p_T invisible 

KEY: TH1F hist10622;1 Mis: pT miss (true - reco)/reco 

KEY: TH1F histl0711;l BJet: b-jets multiplicity 

KEY: TH1F histl0721;l BJet: b-quarks HP multiplicity 

KEY: TH1F hist10723;1 BJet: delta r bjet-bquark HP 

KEY: TH1F hist10724;1 BJet: pTbjet/pTbquark HP 

KEY: TH1F histl0811;l CJet: c-jets multiplicity 

KEY: TH1F histl0821;l CJet: c-quarks HP multiplicity 

KEY: TH1F hist10823;1 CJet: delta r cjet-cquark HP 

KEY: TH1F hist10824;1 CJet: pTcjet/pTcquark HP 

KEY: TH1F histl0911;l Tau: tau-jets multiplicity 

KEY: TH1F histl0921;l Tau: tau-had multiplicity 

KEY: TH1F hist10923;1 Tau: delta r tau-jet, tau-had 

KEY: TH1F histl0924;l Tau: pTtaujet/pTtau-had 

KEY: TH1F histll001;l Calibration: calibration correction 

KEY: TH1F histllOlljl Calibration: pT jets before calibration 

KEY: TH1F histll012;l Calibration: pT jets after calibration 

KEY: TTree ACDTree;! ACDTree 
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This figure "acerdet_structure.png" is available in "png" format from: 


http://arxiv.org/ps/1507.00995vl 



